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X-ray free electron lasers (XFELs) generate sequences of ultra-short, spatially 
coherent pulses of x-ray radiation. We propose the diffraction focusing spectrom- 
eter (DFS), which is able to measure the whole energy spectrum of the radiation 
of a single XFEL pulse with an energy resolution of AE/E «2x 10~ 6 . This is 
much better than for most modern x-ray spectrometers. Such resolution allows 
one to resolve the fine spectral structure of the XFEL pulse. The effect of diffrac- 
tion focusing occurs in a single crystal plate due to dynamical scattering, and is 
similar to focusing in a Pendry lens made from the metamaterial with a negative 
refraction index. Such a spectrometer is easier to operate than those based on bent 
crystals. We show that the DFS can be used in a wide energy range from 5 keV to 
20 keV. 
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Hard x-ray free electron lasers (XFELs) (Emma et al. (2010), 
Ishikawa et al. (2012), Altarelli et al. (2006)) generate 
sequences of ultra-short, spatially coherent pulses of x-ray radi- 
ation. These pulses can be less than 100 fs in duration, and 
have an energy spectrum of the relative width Aw/w ~ 10 -3 
( see for example Altarelli et al. (2006)). Such characteristics 
make XFELs useful for a wide range of applications, including 
potentially 3D imaging of individual macromolecules, looking 
at dynamics of chemical reactions and many others. It is well 
known, that the XFEL pulse energy spectrum has a fine struc- 
ture, which changes chaotically from pulse to pulse. This fine 
structure has a characteristic width on the order of Alu/uj ~ 
10 -6 for hard x-ray energies. Therefore, it is of interest to elab- 
orate a spectrometer with sufficient resolution to resolve these 
fine features for each pulse individually. 

Since each pulse of the XFEL has a unique structure of the 
energy spectrum, the spectrometer has to be able to detect the 
whole spectrum of a single pulse simultaneously. Obviously, the 
scan of the energy region by the monochromator, as it is conven- 
tionally performed at synchrotron sources, can not be applied 
at XFELs. Most of the known spectrometers (for example, 
Kleimenov et. al. (2009), Dickinson et. al. (2008)) of such type 
have the maximum energy resolution of 10~ 5 . They are based 
on the Bragg diffraction of the curved or flat single crystals, 
and their resolution is limited by the angular divergence of the 
reflected beam due to the dynamical scattering effects. Recently, 
a spectrometer based on a curved crystal in the Bragg scatter- 
ing geometry for applications at XFEL sources was developed 
(Zhu, 2012). It still has a resolving power of 2.4 x 10~ 5 . 

We propose a new spectrometer based on the effect of a single 
crystal diffraction focusing (Afanasiev & Kohn, 1977). A sketch 
of the diffraction focusing spectrometer (DFS) is shown in Fig. 
1 . In this spectrometer the incoming divergent beam from a sec- 



ondary source is focused by a single crystal plate at each energy 
near the Bragg angle. Since the Bragg angle depends on the 
energy, x-rays with different energy will be focused at different 
points just behind the crystal under conditions of focusing. Due 
to this property a higher resolution spectrometer can be devel- 
oped. 



Detector 



from 

Source v—K-L 




Figure 1 

Schematic view of the beamline with the diffraction focusing spectrometer. The 
incoming beam from the source is focused by the compound refractive lenses 
(CRLs) and creates the secondary source. The flat crystal is positioned in a Laue 
geometry behind this secondary source. The high resolution detector is located 
close to the crystal in the direction of the diffracted beam. 



In order to achieve high resolution, the effective size of the 
secondary source should be decreased to about few microns 
and the angular divergence of the beam increased to the neces- 
sary value. This can be done, for example, by positioning com- 
pound refractive lenses (CRLs) (Snigirev et al., 1996; 1998) 
upstream from the DFS. Planar Silicon CRLs allow one to 
achieve small focusing distance (few centimeters), and obtain 
the source image of the small size even at high photon energies. 
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An effect of single crystal diffraction focusing was first pre- 
dicted by Afanasiev & Kohn (1977) and considered in detail 
by Kohn (1979). It was experimentally confirmed by Aristov et 
al. (1978, 1980, 1982), and by Koz'mik & Mikhailyuk (1978). 
Recently the effect of focusing in bent-crystal systems was con- 
sidered theoretically by Mocella et al. (2004, 2008), and by 
Nesterets & Wilkins (2008) as a way for improving polychro- 
matic focusing. 

The effect of diffraction focusing is similar to one predicted 
in metamaterials with a negative refractive index for a long 
wavelength radiation. For these materials Pendry (2000) has 
introduced the lens in the form of a plane slab made from 
the metamaterial. The divergent incident beam is focused, first 
inside the slab, and then behind it. 

For x-rays an effective negative refraction occurs for a wave, 
for which, due to dynamical diffraction in the Laue case, the 
Pointing vector has an opposite direction to the angular devia- 
tion of the incident wave from the Bragg angle. As a result, a 
slightly divergent beam incident at the Bragg angle at the single 
crystal is focused, first inside the crystal, and then once again 
behind the crystal. The difference from the Pendry lens is that 
the angular width of the focused beam is equal to the angular 
range of the dynamical diffraction (tens of microradians), and 
is rather small. Nevertheless, this effect allows one to obtain 
a beam of few microns width. The width of the focused beam 
determines the resolution of the DFS. In this work we investi- 
gate the optimal conditions for such a spectrometer based on 
diffraction focusing and show that the energy resolution of DFS 
can reach 10~ 6 . Such resolution will be sufficient to resolve the 
fine structure of the hard x-ray FEL pulse spectrum. 



2. Theory and Analytical Estimates of the Achievable 
Resolution 

Our goal is to obtain an estimate of the resolution of the spec- 
trometer based on diffraction focusing. To get this value it is 
necessary to calculate the full width at half maximum (FWHM) 
of the intensity profile of the focused monochromatic beam at 
the detector position. Main components of the DFS are shown 
in Fig.2. The crystal in the form of plane parallel plate of thick- 
ness t is located at the distance Lq from the point source located 
at the coordinate x s in the plane of the secondary source formed 
by the CRL (Fig. 1). The detector plane is positioned at the dis- 
tance L\ and is perpendicular to the direction of the diffracted 
beam with the wave vector kh. The coordinate xj is related to 
the point of observation at the detector. The angles fo and 
between the normal to the crystal surface no and the wave vec- 
tors ko, kh determine the asymmetry of the diffraction; their 
sum is equal to 29 B , where 9 B is the Bragg angle. 

The amplitude of the diffracted wave Eh{xd) at the detector 
position in the scheme shown in Fig. 2 was obtained by Afana- 
sev & Kohn (1971). It was derived by introducing the crys- 
tal propagator in the form of the Green function for the Tak- 
agi equations (Authier, 2005), and the Kirchhoff propagator, 
describing the transmission in free space in paraxial approxi- 
mation. In the following paper by Kohn et al. (2000) the final 
form of this amplitude was obtained by using Fourier transform 
in reciprocal space. Using this approach the amplitude of the 
diffracted wave Eh{xd) can be presented in the following way 

oo 

E h (x d )=C 1 J ^±J2±[±^M^±(V)-M ± (v))}. (1) 
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Figure 2 

Main components of the diffraction focusing spectrometer. The secondary 
source plane is described by the coordinate x s , crystal is positioned at the dis- 
tance Lq from the source, detector is located at the distance L\ from the crystal 
in the direction of the Laue diffracted beam defined by the vector kh. The detec- 
tor plane is described by the coordinate xj. 
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In these equations the following parameters are used: A is the 
wavelength, is the absorption coefficient, L, = Lq + L\, 
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In Eqs. (2-3) A is the extinction length, fo is an effective 
focusing thickness of the crystal, L is an effective distance from 
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the source to the detector, xo and xi are the displacements 
of the image due to asymmetry of the diffraction geometry, 
Xo,h — XrOM + iXiO,h are the Fourier components of the complex 
susceptibility of the crystal for the angles of scattering and 
28b- The intensity distribution at the detector position Xj for a 
point source at the position x s is given by the square modulus of 
the amplitude Eh(xd) (Eq. (1)), Ih{xd) = \Eh(xd)\ 2 . For geomet- 
rical optics description of the effect of diffraction focusing see 
Appendix. 

The DFS energy resolution is determined by the condition 
that a shift of the photon energy AE leads to a shift of the 
diffracted beam center by x^ at the detector plane. It was 
shown by Kohn et al. (2000) that 
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We will consider now the case of the crystal thickness t = to, 
when focusing conditions are satisfied at the detector position. 
To obtain an estimate of the DFS resolution the shift of the beam 
position in equation (4) has to be compared with the beam 
size (FWHM) A A at the detector position. To obtain an analytical 
estimate of the beam size we will consider two cases of weakly 
absorbing (pok) <C 1) and strongly absorbing (poto 3> 1) crys- 
tals. 

In the first case, we can neglect absorption and consider 
in equation (1) the central part of the integration region with 
77 -c 1 ■ This approximation is justified if the additional condi- 
tion A <C to is satisfied as well. Then, we can consider only one 
term in the equation (1) with the positive sign, expand S(rf) in 
the Taylor series up to the fourth order in T) and obtain 



Eh(xd) ~ C\ I dr/exp inrjs + i 
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is the dimensionless coordinate at the detector position. As soon 
as we consider small values of 77, the integral in the equation (5) 
reaches its maximum value if the following conditions t = to, 
and s — Q are satisfied simultaneously. In fact, the first condition 
is the condition of focusing. It also defines the focusing crystal 
thickness fo, that is according to its definition in equation (3), 
is proportional to the effective distance from the source to the 
detector L. This effective distance is equal to the total distance 
L t in the symmetrical (fi = 1) scattering geometry. 

The integral in equation (5) is similar to the one describing 
the diffraction by a parabolic mirror (Berry & Klein, 1996). 
Therefore, the diffraction pattern in the case of small absorp- 
tion, when the phase term plays the main role, is similar to the 
one resulting from the parabolic mirror focusing. 

To obtain an estimate of the focused beam size, we note that 
the exponential in the integrand becomes a fast oscillating func- 
tion for small values of s outside the region where the term rf is 
less than tt. Therefore, we can exclude this region from the inte- 
gration, and consider the integral in the finite limits from —770 



to 770, where 
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Thus, we arrive to the conclusion that the amplitude of the 
focused wave in a weak absorption case is described by the 
sinc-function, for which the size (FWHM) of the intensity dis- 
tribution is given by 
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Figure 3 

The case of a Silicon 220 reflection, (a) Two-dimensional distribution of the 
intensity, (b) Intensity in the central position of the image xj = as a function 
of the crystal thickness t. 

The second case of strong absorption was discussed earlier by 
Afanasev & Kohn (1977). Again, in equation (1) we take into 
consideration only one term with the positive sign for which the 
Borrmann effect takes place in the region of small 77 <C 1. In 
this case the absorption term can not be neglected, and we can 
expand functions <t> + (77) and M + (77) in the exponent of the Tay- 
lor series up to the second order in 77 and obtain for the ampli- 
tude 
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The integral (10) is similar to the well known Fourier trans- 
form of the Kirchhoff propagator 

00 2 

/ t exp { iqx - V 2 ) = WF exp (12) 

—00 

and can be calculated analytically as 

Eh(x d ) « C 3 exp (i2-K(T2S - 1 

C 3 = ^(^"exp^^-). (13) 

Now the size (FWHM) of the intensity distribution in the focus 
for a strong absorption case is equal to 
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Direct computer simulations show that both values of the 
peak size A xp and A xa , defined by the equations (9) and (14), 
give an excellent estimate within the region of their validity. 
However, in some intermediate cases the region of integration is 
limited by the phase but absorption can not be neglected inside 
this region. In these cases there is no analytical solution for the 
amplitude of the diffracted wave and for the estimate of the 
beam size one should calculate the following integral numer- 
ically 



i/ii 

Ehixd) = C2 / dr]exp(msr) - cnrj 2 + azrf), (15) 



and then estimate A x as the FWHM of the intensity curve. 

Finally, the energy resolution of DFS can be obtained using 
the Rayleigh criterium (Born & Wolf, 2000), i.e. two peaks are 
resolved, if the intermediate point between them has a height 
less than 80% of the maximum height. Taking into account 
equation (4) we obtain for the energy resolution 



AE 
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(16) 



To increase resolution one has to consider Lo to be as large as 
possible, whereas L\ considered to be as small as possible. Ana- 
lytical results obtained for the beam size A x (equations (9) and 
(14)) allow one to estimate the influence of the asymmetry fac- 
tor p. 



Figure 4 

The case of a Diamond 220 reflection (a) Two-dimensional intensity distribu- 
tion, (b) Intensity in the central position at xj = as a function of the crystal 
thickness t. 

In the case of low absorption substituting equation (3) into 
equation (9) and considering the distance L\ = we obtain for 
the beam size (FWHM) 



_ 0.83 A J/ X 2 (2L () ) 

ixP ~ 27T \ Xrh \ m p m 



\l/4 



(17) 



Substituting this equation in the expression (16) we obtain that 
the resolution can be, in principle, increased in the asymmet- 
ric case with ft <C 1. Unfortunately, the gain is very small since 
AE/E oc f3 m . It is also known that small values of (3 are difficult 
to realize in practice because the surface of the crystal must be 
cut at the definite angle with the reflecting atomic planes close 
to the Bragg angle, which also depends on the radiation energy. 
On the other hand, as it follows from Eqs. (2) for the small val- 
ues of j3 the source of the size S will increase the size of the 
focused beam to the value Sift. This will reduce an effective 
resolving power of the DFS. In the following, we will perform 
computer simulations and present results only for the symmet- 
rical (fi = 1) scattering geometry. 

3. Results and Discussion 

In practical applications of diffraction focusing one should have 
single crystals of sufficient thickness. Since large single crystals 
of the highest quality are obtained from Si and Ge we consid- 
ered these two cases in our simulations. Taking into account 
the high power of the XFEL pulses and a demand of high ther- 
mal conductivity, we also considered a single Diamond crystal, 
below referenced as C, in our simulations. Only strong 220 and 
400 reflections were considered. The detailed computer simu- 
lations were performed for the photon energy E — 12 A keV 
(A = 1 A). At this photon energy, the absorption coefficients 
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are: /in = 39.8 cm -1 for Si, /in = 810.0 cm~' for Ge, and 
/io = 3.15 cm" 1 for C. We note that at this photon energy and 
with crystal thicknesses defined by the parameter to in equation 
(3) Si and Diamond crystals can be considered as thin ones, 
whereas Ge crystal is of the intermediate thickness. 
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Figure 5 

Comparison of the intensity profiles at the detector position, (a) Silicon 220 
reflection, crystal thickness to = 360 fim. Direct calculations using Eq. (1) 
(red line) and analytical results based on Eqs. (8) (blue line) and (15) (green 
line) are presented, (b) Diamond 220 reflection, crystal thickness tg = 158 fj,m. 
Direct calculation using Eq. (1) (red line) and analytical results based on Eqs. 
(8) and (15) (blue line) are presented. Both intensity profiles obtained from the 
analytical result completely coincide with each other in this case. 

It follows from equation (16) that for a best DFS resolution, 
the distance Ln should be as large as possible while the distance 
Li, on the contrary, should be as small as possible. In our sim- 
ulations we considered the following distances Lq = 10 m and 
L\ = 0.1 m. Our simulations show that small variations of these 
values does not significantly influence results. 

The intensity distribution at the position of the detector as a 
function of the crystal thickness t for the Laue diffracted beam, 
Si 220 reflection, and the point source is presented in Fig. 3. 
The two-dimensional (2D) intensity distribution in the plane of 
parameters (x^, t) is presented in Fig. 3(a) and a cut through this 
distribution at the center of the detector position Xd — in Fig. 
3(b). Intensity values are normalized to the values at the same 
position along the optical axis but without a crystal. It is well 
seen in this Figure that at the crystal thickness t = to = 360 
/im the incoming x-ray beam is focused on the exit surface of 
the crystal. The relative intensity is also increased to about three 
times comparing to a beam without a crystal. What is also well 
seen in Figure 3(a) are oscillations of intensity as a function 
of crystal thickness. They were observed before in experiments 
on the wedge shaped crystals (Aristov et al., 1978; 1980). For 
the small thicknesses t < to, when absorption can be neglected, 



these oscillations are due to a strong interference between two 
types of Bloch waves (corresponding to two signs in equation 
(1)) inside a crystal. The form of these oscillations is different 
from the ones observed with the point source at the entrance 
surface of the crystal and described by Kato (1961), this is why 
they are called anomalous Pendellosung effect. 

The period of oscillations at xj =0 (see Figure 3(a)) is 
equal to the extinction length A (see eq.(3)). In our scattering 
conditions A = 25.1 /im. On the other hand, the thickness of 
diffraction focusing to = 360 fim, so the condition A <C to 
is well satisfied. With increasing crystal thickness the weakly 
absorbing field becomes larger due to focusing effect whereas 
the strongly absorbing field becomes smaller due to both defo- 
cusing and absorption effects. However, for the focusing crystal 
thickness the interference between two fields still exists, and it 
leads to oscillations of the peak intensity. For the conditions of 
the optimal resolution the crystal thickness in the focal region 
should be considered that provides the minimum beam size. 

The results of a similar computer simulation for the case of 
Diamond 220 reflection are presented in Fig. 4. In this case 
the crystal thickness at which focusing occurs to = 158 fim 
is smaller than in a Si crystal, while the extinction length is 
slightly larger A = 30 /im. The absorption is also smaller 
for a Diamond crystal. As a result, strong interference fringes 
do not vanish at the crystal thickness corresponding to focus- 
ing conditions. Still, the peak with the smallest width can also 
be observed in this case. It is interesting that in this case for 
the large crystal thicknesses well pronounced Kato oscillations 
can be observed. They have the opposite sign of curvature 
comparing to the anomalous Pendellosung fringes discussed 
before. One can call fringes at small crystal thicknesses as plane 
wave interference, and fringes at the large crystal thicknesses as 
spherical wave interference. The latter occurs at crystal thick- 
nesses just behind the crystal thickness fn corresponding to 
focusing conditions, where the focus itself can be considered 
as a secondary source. It is similar to the case considered by 
Kato with the source at the entrance surface of the crystal, if 
we separate the crystal into two parts. The first part creates the 
secondary source and the second part gives the spherical wave 
interference. Our simulations show that the Diamond crystal 
can be considered as a good candidate for a DFS because, 
despite the interference effects, one can observe the focused 
beam with the smaller size. 

We also performed simulations for Si 400 and Diamond 400 
reflections and obtained the crystal thicknesses to correspond- 
ing to focusing conditions 182 /im and 80 /im, respectively. 
However, in these cases we obtained even stronger interfer- 
ence effects. For practical reasons it is better to use crystals 
and reflections where these interference effects are minimized, 
due to the fact that oscillations of intensity in the focal region 
smear the focal point and make application of these crystals 
and reflections for DFS more problematic. We also considered 
the cases of Ge 220 and Ge 400 reflections where absorption 
is much stronger than in Si and Diamond crystals, and inter- 
ference fringes are not observed. Unfortunately, the intensity 
of the Laue diffracted beam is also strongly suppressed for Ge 



research papers 



crystals at these reflections making Ge crystals inappropriate 
for DFS applications. Finally, our analysis has shown that for 
DFS applications the most prominent candidates are Si 220 and 
Diamond 220 reflections. 
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Figure 6 

Resolution of the DFS as a function of the photon energy for the case of Si 220 
reflection. Simulations were performed for the crystal thickness t = 360 fim. 

As it follows from our simulations (see Figs. 3(a) and 4(a)) 
the size of the beam in the focal position is of the order of few 
microns. The beam profile in the focal position (at crystal thick- 
ness t = to) for two cases of Si and Diamond crystals (220 
reflection) is shown in Fig. 5. In this Figure results of computer 
simulations are presented together with the results of analytical 
estimates obtained in the previous section. In the case of a Si 
crystal (Fig. 5, (a)), analytical results obtained by equations (8) 
and (15) are slightly different but give an excellent estimate of 
the beam shape in the focus comparing to the curves obtained 
from the computer simulations. In the case of a Diamond crys- 
tal (see Fig. 5 (b)), both analytical curves completely coincide 
with each other in the focal region. It is not surprising as both 
estimates should give good results for a low absorption, that 
is the case for the Diamond crystal. The approximate curve is 
close to the one calculated according to the expression (1), how- 
ever, it has slightly larger width. In Table 1 the size of the beam 
(FWHM) in the focal position for a crystal thickness t = to, 
obtained from computer simulations in the case of Si and Dia- 
mond crystals, is presented. As it follows from our simulations 
typical beam sizes are on the order of few microns, that means 
that the detector should have resolution on the order of a micron, 
or better. 



Table 1 . FWHM of the beam in the focal position, (/im) 





Si 220 


Si 400 


C220 


C400 


Eq. (1) 


6.6 


7.9 


8.2 


10.7 


Eq. (8) 


6.7 


8.9 


9.6 


12.7 


Eq. (15) 


7.1 


9.3 


9.7 


12.8 



The energy resolution of the spectrometer for the photon 
energy of 12.4 keV, reflections 220 and 400 for Si and Dia- 
mond crystals was estimated using Eq. (16) and is summarized 
in Table 2. Our results show that at this energy the resolution 
of DFS is within the range (2 - 3) x 10~ 6 . As it follows from 



our simulations, the highest resolution is achieved with the Dia- 
mond crystals and 400 reflection. It is interesting to investigate 
energy resolution as a function of energy at the fixed crystal 
thickness. Results of our simulations for the energy range of 
5 keV to 20 keV, Si crystal, 220 reflection, and crystal thick- 
ness t = 360 /im are presented in Fig. 6. As can be well seen 
from this figure, at lower energies about 6 keV resolution is 
even higher than at the photon energy 12.4 keV. At higher ener- 
gies resolution decreases with the energy and starts to oscillate. 
However, up to the energy of 20 keV it stays below the value 
5 x 10~ 6 that is sufficient for the effective use of DFS. The pres- 
ence of oscillations in the energy dependence of the resolution 
can be understood due to the fact that the focusing crystal thick- 
ness to is energy dependent. As a result, interference between 
two fields (strongly and weakly absorbing) in the crystal can be 
observed in the energy range as well. 



Table 2. The energy resolution of DFS 
at 12.4 keV, (x 10 6 ) 



Reflection 


Si 


C 


220 


3.1 


2.5 


400 


2.7 


2.0 



4. Applications of DFS at XFEL 

Here we will discuss in more detail whether the resolution of the 
DFS described in the previous sections is sufficient to resolve 
the fine structure of individual XFEL pulses. To answer this 
question we simulate the spectral distribution of the incoming 
XFEL pulse and analyze its shape after Laue diffraction from 
the crystal. To simulate the initial XFEL energy spectrum we 
used an approach proposed by Pfeifer et al. (2010) when a dou- 
ble Fourier transform in the spectral and time domains is used 
starting with the spectrum containing 'white' noise in the phase. 
The initial XFEL pulse generated according to this procedure 
is shown in Figure 8(a). We assumed here an incoming pho- 
ton energy of 12.4 keV, pulse duration T = 100 fs, and spectral 
width Aw/w = 10~ 3 . As it is well seen from this figure, the gen- 
erated pulse spectrum has a typical behavior of a SASE (self- 
amplified spontaneous emission) XFEL pulse (see, for exam- 
ple, Altarelli et al. (2006)). The spectrum of this pulse obtained 
after Laue 220 diffraction from a Si crystal is presented in Fig. 
8(b). In our simulations we considered the experimental geom- 
etry shown in Fig. 2 with the distances Lq = 10 m and L\ =0.1 
m. It is well seen in Fig. 7 that all main features of the energy 
spectrum of the incoming XFEL pulse are well resolved and 
reproduced at the detector. However, due to a finite resolution 
of DFS, a small additional background can be observed in Fig. 
8(b) as well. 

In the simulations presented in Fig. 7 a point source and a 
wide angular divergence of the incident beam was assumed. In 
principle, the finite source size could reduce the energy reso- 
lution of the DFS. As it follows from equation (2), in order to 
keep the resolution value of DFS on the same level, the source 
size should be significantly smaller then the beam size in the 
focal position. 
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Figure 7 

(a) Simulated XFEL spectrum with the incoming photon energy of 12.4 keV, 
pulse duration T = 100 fs, and spectral width Au>/u> = 10~ 3 ; (b) Correspond- 
ing intensity distribution at the detector after 220 diffraction of this pulse from 
a Si crystal. Insets show an enlarged part of the spectrum. 

According to our results obtained in the previous section, the 
source size on the order of, or smaller, than one micron will 
be sufficient to keep resolution on the level of 10 -6 . Since at 
XFELs the source sizes are on the order of few tenth of microns 
(Altarelli et al. (2006)), the focusing optics is an important ele- 
ment for an effective use of DFS. 

Another important point for an effective use of DFS at XFELs 
is the possibility to measure the whole energy spectrum of a 
single pulse. Two conditions should be satisfied in order to 
reach this goal. The first one puts limitations on the angular 
divergence 9 of the incoming beam. Taking into account that 
divergence of the incoming beam should be bigger than the 
acceptance angle defined by the Bragg condition we obtain 
9 > tan 9b(Aoj/uj). For a typical bandwidth of XFEL sources 
Au)/lu ~ 10~ 3 and reflections considered in our simulations 
with tan 9b ~ 0.3 -j- 0.7 we obtain the following condition for 
a minimum divergence 9q > 7 x 10~ 4 . The angular aperture 
9/ of the parabolic absorbing lens is equal to 9/ = S(X/R/3) m 
(Kohn, 2012), where the complex index of refraction is assumed 
as n — 1 — S + ij3, and R is the curvature radius at the apex of 
the parabola. For Si lens and A = 1 A absorption and refraction 
indices are equal to S = 3.2 x 10~ 6 , /3 = 3.2 x 10~ 8 . Therefore 
from the inequality 9/ > 9 we arrive to the condition R < 0.8 
fim that is easy to satisfy in practice. 

The second condition follows from the fact that the detec- 
tor window and the crystal longitudinal size have to be larger 
than the spectrum size Axj. It follows from Eq.(4) that Ax</ > 
LqOq = 7 mm for Lrj = 10 m and 9q = 7 x 10~ 4 . For a pixel 
size 5 fim the detector has to have more than 1400 pixels. This 
can be achieved in modern detectors. 



5. Summary 

In summary, we propose a spectrometer for XFEL based on 
the effect of diffraction focusing, that theoretically allows to 
achieve energy resolution about AE/E ~ 2 x 10~ 6 . This is better 
than in most existing spectrometers, and should allow, in prin- 
ciple, to resolve the fine structure of individual XFEL pulses. 
We investigated different crystals and possible reflections and 
came to the conclusion that the favorable ones are Silicon and 
Diamond single crystals with the 220 reflection. We obtained 
analytical solutions for the spectrometer resolution in two lim- 
its of a strongly and weakly absorbing crystals and compared 
them with the computer simulations. In order to measure the full 
spectrum the incoming beam should have sufficient divergence 
that can be obtained, for example, by focusing of the incoming 
beam with CRLs. The proposed spectrometer is easy to oper- 
ate and can be installed at any hard x-ray beamline at XFEL. 
By analyzing the resolving power of the DFS we determined 
that detector resolution should be on the order of one micron, 
or below what can be achieved by present detectors. Finally, we 
foresee that spectrometer based on diffraction focusing will be 
widely used for the spectral analysis of individual pulses of the 
present and future hard x-ray FELs. 
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Appendix A 

Geometrical Optics Approach for Description of 
the Diffraction Focusing 

Here, we will describe a process of diffraction focusing on 
the basis of the geometrical optics approximation (see Fig. 8). 
In geometrical optics the ray trajectories from the source to the 
detector are considered. First, the base trajectory, which is used 
in the derivation of the equation (1) has to be defined. This 
trajectory starts from the source and goes to the crystal at the 
Bragg angle with the diffracting atomic planes, after taking into 
account the refraction at the crystal surface. Inside the crys- 
tal it goes along the atomic planes in the symmetrical case of 
diffraction. Behind the crystal it goes along the direction of the 
reflected beam. This trajectory is shown in Fig. 8 by the red 
line. It defines the origin on the detector coordinate Xd axis, i.e. 
it corresponds to the point x ( i = 0. The integration variable rj in 
Eq. (1) is proportional to the angular deviation from this base 
trajectory. Each value of 77 corresponds to the ray, and all rays 
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make a total contribution the detector image. As usual, the tra- 
jectory of the ray is defined from the calculation of the integral 
(1) by means of the stationary phase method. 




2d focus 



of the dispersion surface which are well known in the theory of 
plane wave diffraction in crystals (Authier, 2005). 

We are interested in the upper sign, i.e. weakly absorbing 
wave. For this wave we obtain the following trajectory equation 
for small 77 <C 1 



Xd 
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(20) 



As it follows from Eq. (20), the coordinate Xd decreases with 
increasing the crystal thickness t . At a certain crystal thickness 
t = to the coordinate Xd stays zero for all values of 77 inside the 
region where 77 <C 1 . At this thickness all rays go to the same 
point, that means focusing. If the crystal thickness exceeds that 
focusing distance t > to, then focusing takes place inside the 
crystal (see Fig. 9), and rays become again divergent inside a 
crystal, but now xj is negative for a positive 77, and vice versa. 

Let us consider now the case t > to and the distance L\ behind 
the crystal. We obtain 
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(21) 



Here the first term in the square brackets is negative, but the sec- 
ond term is positive. Therefore, the second focusing is possible 
at the definite distance behind the crystal (see Fig. 8). 



Figure 8 

Schematic view of diffraction focusing based on geometrical optics approach. 

The main contribution to the integral corresponds to the 
region where the condition 



d®±(y) 







(18) 



is satisfied. We are interested in the small values of 77 <C 1. 
Therefore, we can consider only the linear term in Taylor series 
over 77. In addition, we restrict ourselves by the symmetrical 

case j3 = 1. 

Let us consider, first, the case in the front of the crystal. Then 
L\ = t = and from Eqs. (18), (2) and (3) we have 



Xd 



\Xrh\ 

S B 



L T). 



(19) 



As it follows from Eq. (19), the coordinate Xd is proportional to 
the angular deviation parameter 77. It is the usual situation for a 
spherical wave. Let us consider, in addition, the crystal plate of 
definite thickness t. Now we obtain two trajectories according 
to the two signs in the equation. They belong to two branches 
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